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SUMMARY 


The influence of transverse shear stresses on the fatigue crack growth rate in thin 2024-T3 
alummum alloy sheets is investigated experimentally. The tests are performed on double-edge 
cracked sheets in cyclic tensile and torsional loading. This loading generates crack tip stress in- 
tensity factors in the same ratio as the values computed for a crack lying along a lap joint in a 
pressurized aircraft fuselage. The relevant fracture mechanics of cracks in thin plates along with 
the details of the geometrically nonlinear finite element analyses used for the test specimen cali- 
bration are developed and discussed. Preliminary fatigue crack growth data correlated using the 

U y ^ pled ® tress intensity factor calibration are presented and compared with fatigue crack 
growth data from pure A Kj fatigue tests. 


INTRODUCTION 

Recently the problem of a crack near a lap joint (see figure 1) in the pressurized fuselage 
of an airplane was recognized as a situation in which material near the crack tip is subjected to 

0 ensile and transverse shear (i.e., out-of-plane tearing) stresses. Concern over the effect of 
earing stresses on the fatigue crack growth rate prompted researchers at NASA Langley to spon- 
sor a program of analytical, numerical and experimental research at Cornell University aimed at 
measuring fatigue crack powth rates in such situations and developing the theory and numerical 
procedures necessary to interpret and to apply the experimental data. 

Pr0 ^ em is J! lustrated in 1- Crack tip tensile stresses arise from the hoop 

s in the fuselage skin, while the out-of-plane tearing stresses arise from the pressure differ- 
ence pushing out on the skm. As shown in figure 1, one side of the crack remains riveted to 
the stringer and to the adjoining skm forming the lap joint. This side of the crack is much stiffer 
than the other side which is only a single sheet thick. This less stiff side bulges out relative to 
tiae Stitt side, resulting in tearing stresses at the crack tip. 

The stresses at the crack tip resulting from this loading can be described as a combination 
of membrane stresses constant through the thickness and transverse shear or tearing stresses 
that have a through the thickness variation. As defined in figure 2, for general loadings of thin 
cracked plates under membrane and out-of-plane loads, two fracture modes with corresponding 
stress intensity factors Kj and K n can be identified with the membrane loads, and two fracture 
modes with stress intensity factors k x and k 2 can be identified with the out-of-plane loads. 

1 This work was performed with support from NASA Langley, under award NAG-1-1311. 
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Ass uming that the crack tip stresses are a superposition of the stresses from the membrane 
and out-of-plane loads, and that the out-of-plane part can be described using Kirchhoff plate the- 
ory, the stresses on a plane ahead of the crack {0 = 0 in the coordinate system of figure 3) are 
[1-3] 

Kj ki 2x 3 
1722 “ IWr + s/Tt h ’ 

__ Kii , Afr ( 1 + V \ 2^3 x 

ai2 ~ V2^ + V^\3 + vJ h ’ 

-k 2 h [\ / 2x 3 V' 

4723 “ 2(2r) 3 / 2 (3 + v) [ V h ) l 

where u is the Poisson ratio and h is the plate thickness. The in-plane stress components, cr\ 2 
and < 722 , each have a constant and linear thickness variation term resulting from membrane and 
bending loads. The out-of-plane stress component, <r 23 , varies parabolically through the thick- 
ness of the plate and has an r -3 / 2 singularity, a manifestation of the limitation of Kirchhoff plate 
theory to satisfy completely the natural boundary conditions occurring at a free crack surface. 
Note the absence of 1/0 r in the k u k 2 terms. This makes the definitions of fci and k 2 correspond 
to those introduced by Sih et al.[3] and is the form adopted by handbooks of stress intensity fac- 
tors^]. For more details and references on crack tip fields in plates loaded out-of-plane see [5,6]. 

Hui and Zehnder[6] argue that crack initiation and growth in thin plates under combined 
membrane and out-of-plane loads can be correlated with the stress intensity factors Kj , Kj j , k\ 
and k 2 (see eq. (1)). For the lap joint crack it appears that only Kj and k 2 are significant. The 
other two stress intensity factors are zero or nearly zero in this case. The relative importance of 
the tension and tearing stresses can be assessed by computing the Ki and k 2 for a crack located 
along the lap joint and comparing the values. Such a computation was performed by Britt [7] and 
is summarized here in figure 4 where the geometry of the problem is shown, and the stress inten- 
sity factor k 2 is plotted against Kj for a variety of crack lengths. For the crack situated along 
the stringer k 2 becomes comparable to Kj. A second scenario in which k 2 is relatively large oc- 
curs for a curved crack in a fuselage as studied by Potyondy et al.[8]. In this problem, the growth 
of a curving crack is simulated; high values of k 2 with respect to Kj were found to occur. 

In both cases, for longer cracks k 2 becomes significant and thus may be an important fac- 
tor in deter minin g fatigue crack growth rate. Whether this is true or not is one of the principal 
subjects of this research. An extensive literature search revealed no data on fatigue fracture un- 
der such loadings, thus we are undertaking a series of fatigue crack growth rate experiments such 
that the Kj and k 2 stress intensity factors are imposed. 

NUMERICAL PROCEDURES 

To correlate crack growth rate with K.j and k 2 , the capability to compute these stress in- 
tensity factors for both the test specimen and for the aircraft fuselage structure is needed. Nu- 
merical analyses are necessary for this computation since the stress intensity factor handbooks 
do not contain k 2 solutions for finite dimensional plates that would be used for experimentation. 
The procedure used here is to compute components of the energy release rate from finite element 
analyses and then calculate Kj and k 2 from the relations between stress intensity factors and 
these components of the energy release rate. 
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The energy released during crack growth is equal to the work done by the tractions acting 
over the area of crack extension. For a linearly elastic plate of thickness h the energy release 
rate, G, for self-similar extension of a through crack lying in the x\, X 3 plane is[9,10] 

J /-A L r+h/2 

° = a1“o 2hAL J 0 J_ h/2 ‘ 7 *( x i’ e = 0)&u i (AL-x l ,6 = „)dx 3 dx 1 , (2) 

where Am* are the components of relative crack tip displacement for a crack extension of A L 
and the repeated index i implies summation over i = 1,2,3. ' 


From this statement of the energy release rate the relationship between G and the stress in- 
tensity factors can be found[6]. The stress and displacement fields from Williams[l,2l are sub- 
stituted into eq. 2, and after some manipulation and grouping of the membrane terms and the 

bending terms, the following relationships between the G components and the stress intensity 
ractors are found J 
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For a general loading the total energy release rate is the sum of the above, 


i.e.. 


G — G\ + G 2 + G$ 4- G 4 + G 5 . 

The numbering of the components of G is arbitrary; the numbering convention chosen corre- 
sponds to the numbering of the degrees of freedom in the finite element program used for the 
stress intensity factor calculations. If the separate terms of the energy release rate in eq. 3 can 
e found from a numerical analysis, then eq. 3 can be used to calculate all four stress intensity 


To validate the method for the calculation of Kj, K H , k, and k 2 , six test analyses of fi- 
nite cracks m infinite plates were conducted. Figure 5 shows these sample problems along with 
their theoretical solutions taken from refs. 3 and 6. Three of these test cases involve only a sin- 
gle mode; thus, the virtual crack extension technique [11-14] is appropriate to compute the single 
stress intensity factor in each of these problems. This technique gives only a total G value as its 
result; consequently, it must be known a prion what single mode is present. The other three test 
cases involve loadings which generate mixed modes at the crack tip. For these cases, the nodal 
release method [15] is necessary to be able to compute the individual G t s and hence the relevant 
stress intensity factors. The nodal release method computes the components of G on a nodal de- 
gree of freedom basis such that each G t component is related to one of the four fracture modes 
shown in figure 2. 


The computations of figure 5 as well as the test specimen calibrations discussed later were 
performed using the STAGS (STructural Analysis of General Shells) code[16]. The plate model 
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had dimensions of 11.0 by 10.0 by 0.090 thick and contained a centrally located straight crack 
of length 1.0 oriented parallel to the 11.0 dimension. No special crack tip elements were used. 
Instead, the STAGS 410 element[17], a four-noded, six degree-of-freedom per node flat plate el- 
ement, was used for the entire mesh. All elements were square inside a near tip region of mesh 
refinement that extends a distance ^ above and below the crack fine and j behind and \ ahead 
of the crack tip (L is the crack length). The standard element size in this region was 

The results of these sample computations are also given in figure 5. For the three single 
mode cases the virtual crack extension results were all within 3.6% difference from the theoretical 
solutions. For the three mixed mode cases the nodal release results were all within 1.3% differ- 
ence from the theoretical solutions. Based on these results the applicability of these two methods 
for the calculation of the four stress intensity factors is established. 

However, for the stress intensity factor calculations needed for the calibration of the test 
specimen both the virtual crack extension technique and the nodal release technique have sig- 
nificant drawbacks. The virtual crack extension technique by its nature does not permit the com- 
putation of the individual components of the energy release rate, rendering it useless for a mixed 
mode situation such as our testing configuration. The nodal release method, although applica- 
ble to any general mixed mode situation, has the disadvantage of requiring two full finite ele- 
ment analyses to be performed per one Gi computation. This would double the already large 
number of required analyses. Consequently, a third method, the modified crack closure integral 
method [18, 19], was used. This method approximates the displacement field over a crack advance 
of A L by using the displacement field behind the crack of length L, thus requiring only a single 
analysis per stress intensity factor calculation. The results already given that relate components 
of the energy release rate on a nodal degree of freedom basis to the four stress intensity factors 
can be applied to the modified crack closure integral method. 

EXPERIMENTAL METHOD 

To determine the fatigue crack growth rate as a function of the two stress intensity factors 
Ki and double-edge cracked plates were tested under constant amplitude cyclic tensile and 
torsional loadings. The test specimens were machined from 2.29mm thick 2024-T3 aluminum 
plates provided by NASA Langley. This is the same batch of aluminum sheets used by Hud- 
son[20] in his pure mode I fatigue crack growth study in 1969. The test specimen dimensions and 
geometry are shown in figure 6. On each side of the test specimen a 6.35mm long, 0.18mm wide 
starter notch was cut using a foil electric discharge machine. The samples were diamond polished 
to a 3/im finish to improve the optical crack length measurement. 

Ideally, the sample loaded into the gripping fixtures is subjected to uniform vertical dis- 
placements (in the y direction; see figure 6) and to a rigid rotation about the y axis at the bot- 
tom boundary. Any deviation from this condition will change the stress distribution, and con- 
sequently the stress intensity factors and rate of crack growth. To prevent any slipping of the 
test specimen during a fatigue test, the gripping fixtures and face plates each had two channels 
cut along their wide direction, into which were inserted 1.52mm diameter steel rods, as shown 
in figure 6. Sandwiching the specimen between fixtures and face plates and tightening the bolts 
pushes the steel rods into the specimen creating a slight indentation which then ensures a main- 
tainable grip line for the duration of the fatigue test. 
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All of the tests were performed using an Instron 1321 tension-torsion servohydraulic test 

— rr ng “ ? nS10, ‘ a “ d ‘ 0rqUe control - The sample, fixtures and part of the testing mu- 
ch, ne are shown in figure 7. The testing conditions were monitored continuously with an oscillo 

Smm l™ Tt * “ mP , Uter ‘ hat <UgitiZed the W ’ stroke ' ‘-que and rotation s£T A 

lo^dTfiHz Jfh each Starter notch ** eyehag with a pure tension 

te^ileW tnl imnSf Y^ and B = 0.7, where R = P min /P m „, P being the axial 
tensile load. Typically 100,000 cycles were required for the fatigue precracking. 

lasting was performed at R = 0.7 for both the tensile and torsional loads. The tensile loads 

noS fr n h aPP TTT ely 2] f t0 45kN ’ and the tor< l ues ran 6 ed from approximately UNm to 
HONm. A thin plate has a relatively low torsional stiffness, thus the specimen experiences large 

na^lT (on th * ° rder ° f 2 °° for a torque of llONm) severely limiting the testing frequency com- 
p d to a pure tension loading test. Above 1Hz the testing machine actuator could not maintain 

lollTt”reUe7ri\ o% ° f the tCStS W6re Performed at 1HZ Where the loads ““' d »n- 

... < ?. r , ack lengths " ere m «asured approximately every 0.5mm of crack growth and were recorded 
with the corresponding number of elapsed cycles as well as the loading conditions. These crack 

n? frtT tW ° * avemng m * crosco P es , one for each crack, mounted t7th?St- 
? . y ertI cal and horizontal microscope movements were measured with digital travel indi- 

s to give an (x,y) location of the crack tip at each measurement. When setting up the test 
the microscopes were rotated to be normal to the plane of the center of the specim^ when tto 
merni orque was applied. This initial rotation was maintained throughout the test. Fluorescent 

the rSon' 1 wTh^tb e “ SUrfaCe With a Sheet ° f white card st °<* was used to illuminate 
±0 05mm W “ h th,s arrangement crack length could be measured to an accuracy of 


TEST SPECIMEN CALIBRATION 

,, A CaUbrat ; 0 . n is needed to relate experimentally measurable quantities, such as the applied 

, • an > ^ rji to t ^ ie stress intensity factors. For the current experiments, P the net 

tension load T, the net torque and, o, the crack length are inputs to a ^hbration Wtfon ^th 

stress intensity factors, K r and k 2 , as output. What makes the calibration for this testing ar- 
rangemont difficult is that K, is not solely a fimction of P and a and hz is not solely a Stan 
of T and a, on the contrary, AT/ and k 2 are both nonlinear functions of P, T and a Perhaps a 

^ °“ g thl l d f regarding the fr^ure mechanics is that the axial stiffness as well 
as the torsional stiffness of the specimen are not just functions of the plate geometry but func- 
tions of the plate geometry and the current loading of the plate. To approach this aspect of the 
experimental work it was deemed early on that a numerical approach would be needed to com- 

PU 'T: i ‘*|, COTre ? t stress . “tonsity factor calibration. The approach thus employed is a fully geo 
metrically nonlinear plate finite element analysis. ® 

The means used to calculate the stress intensity factors from the results of a finite element 

ellZ H i ‘T “ C " bed “ * previo,,s sect * on - Using a half-plate model (see figure 6 for mod- 
the G d ™h.« and b °T da 7 to”* 410118 ) “d the modified crack closure integral method to obtain 
G, values, geometrically nonhnear analyses were performed for a grid of P, T values at five 
.hscrete crack lengths: 6.35mm, 12.7mm, 19.05mm, 25.4mm and swfmm, ch<»en toTpan the 
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experimental results which are for crack lengths from 6.0mm to approximately 32.0mm. The Ki 
calibration “surface” as a function of P, T for three crack lengths (6.35mm, 19.05mm, 31.75mm) 
is given in figure 8. The corresponding k 2 calibration for the same crack lengths is given in figure 

9. 


Although these calibration plots exhibit many interesting features only a few will be dis- 
cussed here. Notice that for higher torque values A"/ becomes less sensitive to increasing axial 
load (see figure 8). This is especially conspicuous for the shortest crack length. Also notice that 
k 2 at higher torque levels decreases as the axial load is increased (see figure 9). This effect is a 
result of the dependence of the torsional stiffness on the axial tension load in that as the axial 
load is increased, the torsional stiffness of the plate also increases. 

Further results from the calibration calculations although not shown in figures 8 and 9 
also show that Ku and k\ are small compared to Kj and k 2 , respectively, though they are 
present [7, 8]. This result is a welcome one, since the crux of this research depends on the near 
tip stress field in the test specimen being essentially the same as the near tip field in the actual 
cracked aircraft fuselage. Simulations of crack growth in fuselage bays show that indeed the Ku 
and k\ components are small compared to Kj and k 2 , respectively. 

MIXED MODE FATIGUE CRACK GROWTH RATE DATA 

The experiments completed to date cover a wide range of cyclic loadings. The stress inten- 
sity factors corresponding to these tests were calculated based on the calibration curves of figures 
8 and 9 and the maximum values at the start of a test and the end are plotted as the dashed 
lines in figure 4. This figure shows that the experiments cover the lower end of the Ki and k 2 
range for a crack lying along a stringer in a lap joint. 

Results from two of the tests are presented here. Analysis of the remaining twenty tests 
completed to date is in progress. Figures 10 and 11 show the envelopes of stress intensity fac- 
tors for test nos. 14 and 18. In both tests the torque range was the same, but the axial load in 
test no. 14 was ten times larger than the axial load range used for test no. 18. In both tests Ki 
is larger than k 2 for most of the test except for the longest crack lengths in test no. 18. 

Fatigue crack growth rates as well as crack length versus elapsed cycles data for each test are 
shown in figure 12. Shown on the growth plots are benchmark data from Hudson and Newman 
for pure mode I fatigue crack growth[20,21]. 

In each test the crack grew intially at a rate close to the pure mode I rate but slowed down 
dramatically as the crack grew. Thus the ^ data fell well below the reference lines for longer 
crack lengths. T his appears to be caused by roughness induced crack closure. Dining an experi- 
ment it is observed (by eye and ear) that the crack faces are in contact as they slide with respect 
to each other. A s canning electron microscope image of the fracture surface (figure 13) shows the 
region where crack contact has occurred. In this figure the horizontal direction is parallel to the 
plate surface; the direction of relative crack face sliding is vertical. The fracture surface is flat- 
tened out and scoring in the sliding direction is clearly seen. 

However, crack growth rate acceleration does occur towards the end of the test data pre- 
sented in figure 12. It is thought that the crack surface abrasion just mentioned might wear away 
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enough of the contacting surface that at some critical point this effect is minimized and the crack 
growth rate once again increases. If this sort of three-phase behavior (i.e., initial pure mode I 
type growth, then growth rate decrease followed by an increase) is common to all such loading 

situations for the testing configuration is still unknown, but should be elucidated after all of the 
test data is analyzed. 


DISCUSSION 

The experimental data indicate that in the regime of stress intensity factors for which tests 
were performed the crack growth rate is initially reduced by the presence of k 2 resulting from 
crack face contact, but that generally this condition does not persist and crack growth rate does 
eventually increase Whether this is true at higher values of Kj and k 2 is unclear at this time, 
ntuitively it would seem logical that by increasing the value of the stresses at the crack tip the 
crack would grow faster. This may be the case at the higher levels of K z and k 2 where the crack 
faces should separate more and crack face contact should not retard crack growth rate. Thus fu- 
ture experiments may concentrate on higher stress intensity factor levels. 


SUMMARY AND CONCLUSIONS 


Although it is still unclear as to what systematic role, if any, the presence of a k 2 crack tip 
fracture mode might have on the overall fatigue crack growth rate of a crack in a stiffened fuse- 
age structure, it is clear that Characterizing the problem as one of simple mode I fatigue crack 
growth resulting from the fuselage hoop stress would greatly underestimate the Kj field at the 
crack tip. Geometrically nonlinear analyses and resultant stress intensity factor calibrations are 
necessary to describe properly the crack tip stress field in a cracked fuselage lap joint as well as 
any test specimen configuration used to examine the fatigue crack growth behavior in such struc- 
tures. Whether or not the presence of the k 2 mode partially drives the fatigue crack growth or if 
it is simply - the additional AT/ supplied by the tearing or bulging aspect of the crack through the 
geometrically nonlinear effects that drives the crack growth is still not resolved. 

However, this investigation has seemingly raised as many questions as it may have answered 
The question of how to account for crack face contact effectively and what impact it may have on 
the local loading or unloading of the crack tip is unresolved. Furthermore, the more basic ques- 
ion of whether or not the same phenomenon occurs when a long crack in a fuselage lap joint 
is flapping is of primary concern as to the relevance of this work to the actual cracked fuselage 
problem. To complete some fatigue tests at stress intensity factor levels equal to those experi- 
enced by an actual cracked fuselage might be required to answer this question definitively. Fi- 

n y, a quicker and more accurate method of computing all stress intensity factors from a finite 
element analysis is required. 
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Figure 1. Crack in an aircraft fuselage lap joint along a stringer causing crack face 
bulging. 



Symmetric membrane loading, Ki Anti-symmetric membrane loading, Kn 



Symmetric bending, k t Anti-symmetric bending and transverse shear, k 2 

Figure 2. Fracture modes for a thin, cracked plate: Ki, Kji, fci, &2- 
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Frame 




Figure 4. Stress intensity factors, Ki and k-z, for cracks at three positions in a fuselage 
bay. Points are the computed data from Britt [7]. Dashed lines are the max- 
imum stress intensity factors from the tension-torsion fatigue tests described 
herein. 


902 





Norm. Theoretical Values 

Loading: 

Tension (VCE)’ 

Bending (VCE) 

Transverse Shear (VCE) 

Tension, Shear & 

Bending (NR) 2 

Bending & Transverse 
Shear (NR) 

Tension & Transverse 
Shear (NR) 


k l k 2 

( 6 . 0 ) ( 3 . 0 ) 


5.954 (0.8%) 

2.891 (3.6%) 

5.965 (0.6%) 

5.997(0.1%) 2.967(1.1%) 

2.963 (1.3%) 


K I K n 

(1.7725) (1.7725) 


1.754 (1.0%) 


1.765 (0.4%) 1.763 (0.5%) 


1.760 (0.7%) 


’VCE... virtual crack extension technique 
2 NR... nodal release technique 


Figure 5. Six test cases for the validation of the stress intensity factor computation 
techniques: finite length crack in an infinite plate. 
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y * v 




(c) 


Figure 6.^ (a) Test specimen, dimensions in nun. (b) Boundary conditions used in finite 
element model, (c) Detail of specimen gripping fixture. 
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50 " 

Torque, Nm 25 Q Q Tension, kN 


Figure 8. Stress intensity factor Kj versus axial load P and torque T for three crack 
lengths: 6.35mm, 19.05mm, 31.75mm. 
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Figure 12. Crack length, a, versus total elapsed cycles, N, and crack growth rate &- 
versus AKj for test no. 14 (left) and test no. 18 (right). “A” and “B” des- 
ignations refer to each of the two cracks in a given specimen. Hudson’sf20l 

pure mode I data is shown by the solid line; Newman ’s[21] pure mode I data 
is shown by the dashed line. 









Figure 13. Scanning electron micrograph of the fracture surface in a region where crack 
contact has occurred. The region in the upper part of the picture shows no 
contact: the typical ductile fracture surface is seen. In the lower region the 
ductile fracture surface has been flattened and abraded by the crack surfaces 
sliding past each other. 



